Hole-burning experiments within solvable glassy models 
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We reproduce the results of non-resonant spectral hole-burning experiments with fully-connected 
(equivalently infinite-dimensional) glassy models that are generalizations of the mode-coupling ap- 
proach to nonequilibrium situations. We show that an ac-field modifies the integrated linear response 
and the correlation function in a way that depends on the amplitude and frequency of the pumping 
field. We study the effect of the waiting and recovery-times and the number of oscillations ap- 
plied. This calculation will help descriminating which results can and which cannot be attributed 
to dynamic heterogeneities in real systems. 

PACS Numbers: 64.70.Pf, 75.10Nr 



One of the most interesting questions in glassy physics 
is whether localized spatial heterogenities are generated 
in supercooled liquids and glasses. 

In most supercooled liquids, the linear response to 
small external perturbations is nonexponential in the 
time-difference r. Within the "heteregeneous scenario", 
the stretching is due to the existence of dynamically dis- 
tinguishable entities in the sample, each of them relaxing 
exponentially with its own characteristic time. A dif- 
ferent interpretation is that the macroscopic response is 
intrinsically nonexponential. In the glass phase, the re- 
laxation is nonstationary and the dependence in r is also 
much slower than exponential. 

The heterogeneous regions, if they exist, are expected 
to be nanoscopic. The development of experimental tech- 
niques capable of giving evidence for the existence of such 
distinguishable spatial regions has been a challenge for 
experimentalists. 

With non-resonant spectral hole-burning (NSHB) 
techniques one expects to probe, selectively, the micro- 
scopic responses. |Sj The method is based on a wait, 
pump, recovery and probe scheme depicted in Fig. |l|. The 
amplitude of the ac perturbation is sufficiently large to 
pump energy in the sample, modifying the response as 
a linear function of the absorbed energy. The step-like 
perturbation d is very weak and serves as a probe to mea- 
sure the integrated linear response of the full system. The 
large ac and small dc fields can be magnetic, electric, or 
any other perturbation relevant for the sample studied. 
The idea behind the method is that the comparison of 
the modified (perturbed by the oscillation) and unmodi- 
fied (unperturbed) integrated responses yield information 
about the microscopic structure of the sample. On the 
one hand, a spatially homogeneous sample will absorb 
energy uniformly and its modified integrated response 
is expected to be a simple translation towards shorter 
time-differences r of the unmodified one. On the other 
hand, in a heterogeneous sample, the degrees of freedom 
that respond near the pump frequency f2 are expected to 



absorb an important amount of energy and a maximum 
difference in the relaxation (equivalently, a spectral hole) 
is expected to generate around t ~ 1/Q. 
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FIG. 1. Wait, pump, recovery and probe scheme. 

The NSHB technique has been first applied to the 
study of supercooled liquids. The polarization response 
of dielectric samples, glycerol and propylene carbonate, 
was measured after being modified by an ac electric field. 
M More recently, ion-conducting glasses like CKN || , re- 
laxor ferroelectrics (90PMN-10PT ceramics 

) H and spin- 
glasses (5% Au:Fe) || were studied with similar meth- 
ods. The results have been interpreted as evidence for the 
existence of spatial heterogeneities. We show here that 
their main features can be reproduced by a system with 
no spatial structure. We use one model, out of a family, 
that captures many of the experimentally observed fea- 
tures of super-cooled liquids and glasses as, for instance, 
a two-step equilibrium relaxation close and above T c [[|, 
aging effects below T c fjj, etc. The model is the p spher- 
ical spin-glass [||, that is intimately related to the 
mode-coupling model |J. It can be interpreted as a sys- 
tem of N fully-connected continuous spins or as a model 
of a particles in an infinite dimensional random environ- 
ment, ficj] In both cases, no reference to a geometry in 
real space nor any identification of spatially distinguish- 
able regions can be made. 
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In the presence of a uniform field, the model is 



lower limit SI ~ 27r/(8000e) ~ 0.1 to the frequencies we 



H A S ] = ^2 Jll --- 4 3» Sil 

ii < . . . Kip 



N 

h y^sj 



(1) 



The interactions ./^...i are quenched independent ran- 
dom variables taken from a Gaussian distribution with 
zero mean and variance [Jf 1 ,,,i]j = J 2 p\/(2N P ^ 1 ). p is 
a parameter and we take p = 3. Hereafter [ ]j repre- 
sents an average over P[J] and J — 1. The continuous 
variables Sj are constrained spherically 53i=i s 1 = N ■ A 



stochastic evolution is given to s, s,(t) 



£,-(t) with a white noise with = and (£i(f)£i(t')) = 
2TS(t — t'). When TV — » oo, standard techniques lead to a 
set of coupled integro-differential equations for the auto- 
correlation NC{t,t') = J2f=i[( s i( t ) s j{ t '))}j and the lin " 
ear response R(t,t') = Yli=i 3[( s i(t))] j/3$i{t') i with 

8—0 

5i(t') an infinitesimal perturbation modifying the energy 
at time t' according to H — > H — ^\ ^Sj. The dynamic 
equations read ITTJ 



ftC(M') = -2(i)C(M') + ? / dt"C p - 1 (t,t")R(t , ,t / ') 
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p(p - 1) 



dt"C p - 2 (t,t")R(t,t")C(t",t') 



+2TR(t', t) + h{t) \ dt"h(t")R(t', t") , 



(2) 



d t R{t,t') = -z(t)R(t,t') 



dt"C p - 2 (t, t")R{t, t")R{t", t') , (3) 



The Lagrange multiplier z(t) enforces the spherical con- 
straint and an integral equation for it follows from Eq. (^) 
and the condition C(t, t) — 1. In deriving these equa- 
tions, a random initial condition at to = has been used. 
It corresponds to an infinitely fast quench from equilib- 
rium at T = oo to the working temperature T. The 
evolution continues in isothermal conditions. 

In the absence of energy pumping, these models have 
a dynamic phase transition at a (p-dependent) critical 
temperature T c , T c ~ 0.61 for p = 3. When an external 
ac-field is applied, it drives the system out-of-equilibrium 
and stationarity and FDT do not necessarily hold at any 
temperature. The question as to whether the clearcut 
dynamic transition survives under an oscillatory field is 
open and we do not address it here. We simply study 
the dynamics close to the critical temperature in the ab- 
sence of the field by constructing a numerical solution to 
Eqs. (^) and ([|) with a constant grid algorithm of spacing 
e. We present data for small spacings, typically e = 0.02, 
to minimize the numerical errors. Due to the fact that 
Eqs. (0) and (^) include integrals ranging from to = to 
present time t, the algorithm is limited to a maximum 
number of iterations of the order of 8000 that imposes a 



A word of caution concerning the scheme in Fig. [l] and 
the times involved is in order. For the purpose of col- 
lecting the data for each reference unmodified integrated 
response, the sample is prepared at the working temper- 
ature T at to = and let freely evolve during a total 
waiting time t w + 1\ + t r . Depending on T, this inter- 
val may or may not be enough to equilibrate the sample, 
(ti is chosen as t\ — 2irn c /£l with S7 the angular veloc- 
ity of the field that will be used to record the modified 
curve.) A constant infinitesimal probe 6 is applied after 
t w + t\ + t r to measure 

$(r) = f dr' R(t w + ti + t T + r,t w +t 1 + t r + r') . 
Jo 

As an abuse of notation we explicitate only the r de- 
pendence and eliminate the possible t w + t\ + t r de- 
pendence. The modified integrated response <&* is mea- 
sured after waiting t w , applying n c oscillations of du- 
ration t\ = 2nn c /fl, further waiting t r , and only then 
applying the probe 6. The effect of the ac perturbation 
is then quantified by studying the difference: 



(4) 



We have examined A$ at T = 0.8 > T c and T = 
0.59 < T c . We pump one oscillation with hp = 0.1 and 
later check that this field is small enough to provoke a 
spectral modification that is linear in the absorbed energy 
(see Fig. || below). For simplicity, we start by choosing 
t w = t r = 0. In Fig. we show A$ against logr for 
different fl at T = 0.8. All the curves are bell-shaped 
and vanish both at short and long times. In panel a, the 
Sis are larger than a threshold value fi c ~ 1. The height 
of the peak A<E> m = max(A$) decreases with increas- 
ing frequency reaching the limit A$ m = for SI — > oo. 
In addition, the location of the peak t m moves towards 
longer times when SI decreases. In panel b, SI < S7 C and 
the behaviour of the height of the peak is the opposite, 
it decreases when SI decreases and, within numerically 
errors, its position is either independent of SI or it very 
smoothly moves towards shorter times for increasing SI. 
The nonmonotonic behaviour of $ m with SI is a con- 
sequence of the interplay between t Q , the a relaxation 
time, and 27r/Sl the period of the oscillation. The term 

J^''^ dt"h(t")R(t',t") in Eq. (|) controls the effect 
of the field and, clearly, vanishes in the limits SI — » oo 
and r2 — > 0. The inversion then occurs at a frequency Sl c 
that is of the order of 27r/t Q . These results qualitatively 
coincide with the measurements of the electric relaxation 
in CKN at T < T g in Fig 1 a and b of Ref. g. In Fig. | 
we show A$ against logr for different at T = 0.59. For 
all Q we reproduced the situation of panel a in Fig. as 
if SI > Sl c . We have not found a threshold Sl c , that has 
gone below the minimum SI reachable with the algorithm. 
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FIG. 2. Time-difFerence dependence of the distortion A$ 
due to a single oscillation in a log-linear scale. T = 0.8 > T c , 
hp = 0.1, t w = t r = 0. At high pumping frequencies 
Q > Q c ~ 1, shown in panel a, both the height of the peak 
A$ m and its position t m decrease with increasing frequency. 
The dotted (blue), solid (red) and dashed (black) curves cor- 
respond to Q. = 1,2,5 respectively. In panel b, Q < f2 c ~ 1 
and A$ m increases with increasing while t m is almost un- 
changed. The dotted (blue) and solid (red) curves correspond 
to Q = 0.1 and 0.2, respectively. 
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FIG. 3. The same plot as in Fig. | for T = 0.59 < T c . 
For all pumping frequencies explored, the peak moves towars 
shorter times and its height decreases for increasing frequen- 
cies. The dashed (black), solid (red) and dotted (blue) curves 
correspond to f2 = 0.1,0.5 and 1, respectively. 

The maximum modification of the relaxation A$ m in- 
creases quadratically with the square of the amplitude 
of the pumping field h,F, and hence linearly in the ab- 
sorbed energy, as long as Iif < 1. In Fig. || we display 
the relation A$ m oc hp in a log-log scale for the two 
temperatures explored. The amplitude hp = 0.1 used in 
Figs. @ and is in the linear regime. 



FIG. 4. Check of A$ m oc h% in a log-log scale. Squares 
(pink) correspond to T = 0.8 and circles (blue) to T = 0.59. 
In both cases one cycle of an ac-field with Q — 1 was applied. 
The line has a slope equal to 2 and is a guide-to-the-eye. The 
linear relation breaks down beyond 1%f ~ 1. 
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FIG. 5. The normalized maximum distortion for several 
recovery times t r . In panel a T = 0.59 and crosses, diamonds 
and pluses correspond to = 1,2 and 3, respectively. In 
panel b T = 0.8 and crosses, pluses and diamonds correspond 
to Q = 1, 5, 10, respectively. 

The effect of the pump diminishes with increasing re- 
covery time t r . A convenient way of displaying this 
result is to plot the normalized maximum deviation 
A$ m (t r )/A$ m (0) vs flt r . Using several frequencies and 
recovery times, we verified that this scaling holds for 
T = 0.59 but does not hold for T = 0.8, as shown in 
Fig. ^| This simple saling holds very nicely in the re- 
laxor ferroelectric jl and in the spin-glass || but it is 
very different from the Sl-independence of the propylene 
carbonate S. 

Up to now, the effect of a single cycle of different fre- 
quencies has been studied. Another procedure can be en- 
visaged. Since t\ = 2Trn c /Q, we can change ti by apply- 
ing different numbers of cycles n c while keeping fl fixed. 
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In Fig. ^ we show the distortion due to n c = 10, 2, 1 cy- 
cles with £1 = 10 at T = 0.8. The qualitative dependence 
on n c is indeed the same as the dependence on 1/fi: the 
peaks are displaced towards longer times with increasing 
n c (longer tx). This behaviour is similar to the results 
obtained for propylene carbonate in Fig. 11 of Ref. Bb. 
though we do not reach the expected saturation within 
our accessible time window. 



in order to have unequivocal evidence for them. The de- 
tailed comparison of the experimental measurements to 
the behaviour of glassy models with and without space 
will certainly help us refine the experimental techniques. 
Numerical simulations can play an important role in this 
respect. 

10 I . . . ... . . ... i 
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FIG. 6. Effect of several cycles at T = 0.8 for SI = 10. The 
dashed (blue), solid (red) and dotted (pink) curves correspond 
to n c = 10, 2, 1 respectively. 




FIG. 7. Distortion for t w = 8 with dots (pink) compared 
with the one for t w — with full line (red) at T = 0.59. 

Below T c the nonperturbed model never equilibrates 
and the relaxation depends on t w . Indeed, t a is an ap- 
proximately linear function of t w ]7|,|To| and the distortion 
might depend on t w . We compare A$ vs logr for two 
t js in Fig. 0. 

Finally, we checked that the effect of one or many pump 
oscillations on the difference AC = C* (t w +ti+t r +r,t w + 
t\ +t r ) — C(t w + 1\ + t r + t, t w + 1\ + t r ) is very similar to 
the one observed in A$. This observation is interesting 
since it is easier to compute numerically correlations than 
responses. Figure |8| shows the modification observed at 
T = 0.8 and ft > Q c (to be compared to Fig. |2|). 

We conclude by stressing that we do not claim that 
spatial heterogeneities do not exist in real glassy systems. 
We just wish to stress that the ambiguities in the inter- 
pretation of experimental results have to be eliminated 
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FIG. 8. Change in the autocorrelation. T = 0.8 and SI = 1 
(dashed blue), SI — 2 (solid red) and SI = 5 (dotted black). 
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